Scattering Theory of Photon- Assisted Electron Transport 
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The scattering matrix approach to phase-coherent transport is generalized to nonlinear ac- 
transport. In photon-assisted electron transport it is often only the dc-component of the current 
that is of experimental interest. But ac-currents at all frequencies exist independently of whether 
they are measured or not. We present a theory of photon-assisted electron transport which is charge 
and current conserving for all Fourier components of the current. We find that the photo-current 
can be considered as an up- and down-conversion of the harmonic potentials associated with the 
displacement currents. As an example explicit calculations are presented for a resonant double 
barrier coupled to two reservoirs and capacitively coupled to a gate. Two experimental situations 
are considered: in the first case the ac-field is applied via a gate, and in the second case one of the 
contact potentials is modulated. For the first case we show that the relative weight of the conduction 
sidebands varies with the screening properties of the system. In contrast to the non-interacting case 
the relative weights are not determined by Bessel functions. Moreover, interactions can give rise 
to an asymmetry between absorption and emission peaks. In the contact driven case, the theory 
predicts a zero-bias current proportional to the asymmetry of the double barrier. This is in contrast 
to the discussion of Tien and Gordon which, in violation of basic symmetry principles, predicts a 
zero-bias current also for a symmetric double barrier. 



PACS numbers: 73.23.-b, 85.30.Vw, 73.40. Gk, 72.10.-d 
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I. INTRODUCTION 

Photon-assisted tunneling has been of interest since 
the work of Tien and Gordon [Q and Tucker Q . Carrier 
transmission through barriers with oscillating potentials 
has been analyzed to find the traversal time for tunnel- 
ing j|. Recently photon-assisted tunneling has found 
renewed interest in the field of mesoscopic physics stim- 
ulated by theoretical work by Bruder and Schoeller ||] 
and experiments on quantum dots by Kouwenhoven and 
McEuen et al. ||, and by experiments on superlattices 
by the group of Allen et al. ||jD. Typically of inter- 
est P,p|,[7|-|24]] is the zero-frequency current component 
induced in response to an oscillating voltage. Theoreti- 
cal treatments of photon-assisted electron transport of- 
ten assume that the driving field is known and equals 
the external field. However, the long-range Coulomb in- 
teraction will screen the external field and generates an 
internal potential that can be quite different from the 
applied potential. Similarly, since it is the dc-component 
which is measured, one might think that displacement 
currents play no role. However, the dc-component is a 
consequence of nonlinearitics in the conduction process. 
Clearly, in such a conductor, the current has not only 
a dc-component, but also currents at the frequency of 
the oscillating voltage and its higher harmonics. Not 
only are the dc-currents conserved but also the currents 
at the oscillation frequency and at its higher harmonics. 
Consequently a theory is needed in which all frequency 
components of the current are treated self-consistently. 
Such a theory is developed below. It leads to the conclu- 
sion that the photo-current is induced by a non-linear up- 
and down-conversion of the electric-fields (potentials) as- 



sociated with the displacement current. Bessel functions 
are often a hall mark jjj of the discussion of photon- 
assisted tunneling: However, in the self-consistent theory 
discussed below Bessel functions cannot in general de- 
scribe the relative weights of the sideband peaks since the 
argument of the Bessel functions is not invariant under an 
overall potential shift. Since in nonstationary conditions 
charge accumulation occurs and causes induced fields, a 
selfconsistent treatment of the electron-electron interac- 
tions is important. The issues are similar for theories 
and experiments which investigate photon-assisted pro- 
cess not in the dc-current but in its fluctuations p5| , p6| . 
Here we emphasize mainly the average current and ad- 
dress the fluctuation spectra only briefly in an appendix. 

A convenient description of conduction processes in 
mesoscopic systems which incorporates the role of con- 
tacts and permits to investigate directly the phase- 
coherent transmission from one reservoir to another, is 
the scattering matrix approach p?|p8| . The description 
of linear ac-conduction in response to oscillating poten- 
tials and consideration of the long range Coulomb inter- 
action has already been discussed both for the case of 
zero-dimensional systems and for extended sys- 

tems for which one needs to discuss the entire potential 
landscape [ jBlp^l - A review of this subject can be found 
in Ref . 1 3^] . Here we generalize the scattering matrix ap- 
proach to take into account the nonlinear dependence on 
oscillating potentials. First we consider the response of 
the electrons to a potential applied only to the contacts 
of the sample, assuming the internal potential is kept 
fixed. The response to the total potential will, in a sub- 
sequent step, be calculated selfconsistently in random- 
phase-approximation (RPA). The resulting charge and 
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current conserving theory will be used to investigate the 
photo-induced dc-current in a resonant tunneling barrier. 
As function of Fermi energy and frequency we find large 
differences between the induced internal potential and 
the external applied potential, showing that long range 
Coulomb forces are important for photon-assisted tun- 
neling in mesoscopic systems. Furthermore, interactions 
can give rise to an asymmetry between absorption and 
emission peaks, as well as changing the distance between 
peaks from a multiple of photon quanta to a distance 
depending on screening properties. 

Our discussion is complementary to works which model 
interactions based on a Hamiltonian suitable to describe 
Coulomb blockade effects. The work of Bruder and 
Schoeller jl| also considers coupling to a gate and also 
considers displacement currents. In principle, all the 
questions addressed here can be investigated within such 
a framework. The scattering approach used here has the 
advantage that it is not limited to the tunneling regime 
but can also be applied to conductors which are strongly 
coupled to reservoirs (ballistic or metallic diffusive wires, 
etc). The RPA treatment as it is formulated below does 
have the disadvantage that it is not an appropriate de- 
scription in the case when charge quantization effects 
(Coulomb blockade) are important. However, its con- 
ceptual clarity makes the RPA treatment a useful point 
of reference for comparision of different theoretical dis- 
cussions. 

The basic view taken here is the same as that 
used for the discussion of dc-conductances |2i| and ac- 
conductance |2{|. What is needed is the connection be- 
tween currents at the contacts of the structure and the 
voltages at these contacts. Either the currents or the 
voltages can be controlled. As in the discussion of the 
ac-conductance it is necessary to consider not only the 
mesoscopic conductor itself but all nearby metallic bod- 
ies (gates and capacitors) which interact via long range 
Coulomb forces with the mesoscopic conductor. Let a 
label all the relevant contacts. The current at contact 
a can be written in terms of its Fourier components 
I a {nuj). Here n — is the dc component of the cur- 
rents, and n = ±1 are the Fourier components at the 
driving frequency. Nonlinearities lead to higher harmon- 
ics n = ±2,3, ... Similarly, the voltage at contact a has 
the Fourier components V a (nuj). We emphasize that the 
voltage of a contact is only a well defined quantity if lo- 
cal electric fields deep inside the contact vanish. There 
must, therefore, exist a Gauss volume which encloses 
the mesoscopic conductor |^9| . The electric flux through 
this Gauss volume vanishes. As a consequence the total 
charge Q inside the volume is conserved [^J. Charge 
conservation, and current conservation, apply to each 
Fourier component separately. In particular, we must 
have that the total charge within the Gauss volume van- 
ishes at each frequency, 

Q a {nu) = {). (I) 

A theory for which this holds gives currents which depend 



ultimately only on voltage differences. We call such a the- 
ory of electric conductance gauge invariant To be 
definite let us introduce an expansion parameter e. We 
take the Fourier components of the first harmonic V a (u) 
proportional to e and expand the currents in powers of 
e. The second harmonic voltages V a (2uj) describing two- 
photon processes are then proportional to e 2 . Below we 
write the relationship between currents and voltages up 
to second order in e. The expansion coefficients are con- 
ductances <7 a( 3 7 (noj, mu>) which give the current at con- 
tact a in response to a voltage Vp(ntu) at contact (3 at 
a frequency mo and a voltage a contact 7 at a frequency 
muj. The overall dc-current is 

UO) = lt c {{Vp(0)}} + Ig h [0;{V p (0)}], (2) 
i p a h [0;{V p (0)}} = ^gc&[u,-u;{V s (p)}]V(>(u>)v;{w). (3) 

The first term of Eq. (§), /£ C [{V£(0)}], is the direct cur- 
rent that would be measured in the presence of purely 
static voltages Vg(0), (3 — I, 2, .. applied to the different 
contacts of the sample. In the following, for the direct 
current, we retain the full dependence to all orders in the 
static applied voltages. If the dc-current ^Q C [{Vg(0)}] 
is expanded in powers of the applied voltage then the 
terms linear in the applied voltages determine the dc- 
conductance matrix g a p(0) and the terms quadratic in 
the applied voltages are the dc-rectification conductances 
Qa0.y{O), discussed by Christen and one of the authors 
pq , p9| , which determine the leading order non-linearity 
of the dc I-V-characteristic In addition to these 

contributions to the dc-current which characterize the 
purely stationary transport there is now also a contri- 
bution to the dc-current due to the photon-assisted pro- 
cesses, [0; {V/3(0)}]. In particular, to second order in 
the applied ac- voltages Vg(<jj), carriers which emit and 
re-absorb (virtual) photons are determined by the dc- 
photo-conductance g Q( 3 7 [w, — ui\ {Va(0)}] which depends 
in general also on the dc- voltages Vs(0). These conduc- 
tance coefficients represent an up- and down-conversion 
of the first harmonic voltages. 

The current at the frequency of the oscillating poten- 
tial is in general composed both of a particle current and 
of a displacement current. To be brief we call this current 
simply the displacement current. To linear order in our 
expansion parameter it is given by 

/ a (w)=J]^k{V 7 (0)}]^( W ). (4) 



Here expanding g a f3[uj; \V-y(0)}) in the dc voltages yields 
the equilibrium admittance [ pOf of the mesoscopic struc- 
ture g a p(u}) and the dc-ac-rectification conductance 
g a01 (uj;O). 

The current at 2co is 

I a (2w) = £ fla/9 [2w;{V 7 (0)}]^(2w) 
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+ ^g afh [w,w',{Vs(0)}]V f i{u)V r (w) 

Pi 



(5) 



determined by a second harmonic conductance 
g a f3 [2lo ; {Vy(0)}] and a non-linear up-conversion con- 
ductance g a py[u>, to; {Va(0)}] whereby a second harmonic 
current is generated due to a non-linear combination 
of first harmonic voltages. We emphasize that the ex- 
pansion given here can in principle be carried further 
to an arbitrary order in e. Our task is to find explicit 
expressions for the (non-linear) ac-conductances defined 
in Eqs. ((2]- ||). It is useful, to state first a number of 
general properties of these conductances. 

Current conservation holds for each Fourier component 
separately. Furthermore, since we can break off the ex- 
pansion at any order, current conservation restricts each 
type of conductance coefficient in Eqs. An addi- 

tional restriction imposed on these conductance coeffi- 
cients arises due to the fact that a voltage V(nuo) which 
is applied to all contacts simultaneously can not have a 
physical effect. As a consequence the conductances obey 
the sum rules MM 



^g al3 {ku) ='J2g a/3 (juj) = 



(6) 



for k, j 6 JV. Similarly, the second order coefficients obey 
^2g al3l (kuj,juj) = ^2g a/3l (kuj,juj) 

a p 

= ^2 g Q f3y(kuJ, juj) = 0. (7) 

7 

These sum rules guarantee that the final result will de- 
pend on voltage differences only. 




FIG. 1. Conductor connected to two contacts and coupled 
capactively to a gate. 

Eqs. (|| - ||) are completely general and are applica- 
ble to any phase-coherent multi-terminal conductor. We 
now discuss these general relations for the case of a two- 
terminal conductor capacitively coupled to a gate, a sit- 
uation sketched in Fig. [j], in the limit C\ = C2 = 0. 



This simple arrangement permits us already to point to 
the connection between photo-currents and displacement 
currents. We are interested in the photo-current gener- 
ated by a sinusoidal oscillation of the voltage V g (u) at 
the gate. First consider the displacement current. The 
oscillating gate couples with the conductor in a purely ca- 
pacitive manner. Therefore, the g ag [u>; {Vg(0)}] describe 
capacitive currents and we can write g ag [u; {Vg(0)}] = 
—iu>Cag[w, {Vp(Q)}]. We emphasize that this is a global 
transport coefficient which connects the voltage at one 
contact to the current at another contact. As a conse- 
quence, the capacitance coefficients are not of a purely 
geometrical nature but can be strong functions of mag- 
netic field and the dc gate voltage p9| , [3^ , [34U3l| . Thus, 
the current at contact a is determined by 



I a (u) = -iuC ag [u;{Vp(0)}]V g {u). 



(8) 



Next, consider the dc-photo-current generated by this ar- 
rangement 



F a h (O)=g agg [u>-u-,{V (O)}]\V g (u)\' 2 



(9) 



With any of the equations for the displacement current 
we can eliminate the oscillating gate voltage entirely and 
find 



It (0) = 



1 g agg [w,-Lj;{Vf,(p)}} 
^ \C ag [uj;{Vp(0)W 



i a HI 5 



(10) 



Thus to second order in the oscillating voltages, the 
photo-current is directly related to the displacement cur- 
rent. Since the displacement current is not a prop- 
erty of a non-interacting system but is in an essential 
way determined by the long range Coulomb interaction, 
so similarly, the long range Coulomb interaction must 
play an essential role in determining the photo-current. 
Note that the photo-conductance which enters Eq. jlfj| ) 
is also proportional to lo 2 and the photo-current given in 
Eq. (|1C]) therefore has a well defined zero- frequency limit. 

Now we proceed to find explicit expressions for the 
non-linear conductances introduced above. 



II. OSCILLATING CONTACT POTENTIALS: 
EXTERNAL RESPONSE 

We consider a conductor with voltages which oscil- 
late in time applied to the contacts of the sample or 
to nearby capacitors. First we evaluate the response of 
non-interacting particles with the internal potential kept 
fixed. Only the response to the total potential has phys- 
ical meaning, however, these results are needed in the 
next section for treating the problem with interactions. 

The current operator for current incident in contact a 
in a mesoscopic system can be written as [p9[ 



L(t) = dE J dE 1 [ai(E)k a (E') 

- \>UE)b a (E')] (11) 
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where a Q and h a are vectors of operators with compo- 
nents a an and b an . Here a an annihilates an incoming 
carrier in channel n in lead a and b an annihilates an out- 
going carrier in channel n in lead a. Eq. ( |ll"| ) applies 
for frequencies (E — E')/h small compared to the Fermi 
energy. 

The incoming and outgoing waves are related by the 
scattering matrix 39 s a p via, b Q = J^p s a0&p- In a mul- 



tichannel conductor the s-matrix has dimensions N a x Np 
for leads with N a and Na channels. Here, and in the 
following, greek indices run over all contacts of the con- 
ductors. 

Let us now suppose that a potential variation is applied 
to reservoir a. The potential is eU a (t) = eV a (u>) cos ut, 
where V a (oj) is the modulation amplitude. With this 
potential the solution to the single-particle Schrodinger 
equation at energy E in a is 



tp a<n (x, t; E) = 4> a ,n{x\ E)e 



-iEt/h 



£ 



Tiuj 



iliot 



(12) 



where 4> a ,n{x] E) is the wave function describing an in- 
coming (or outgoing) carrier in contact a in channel n 
in the absence of a modulation potential, and J\ is the 
Tth order Bessel function. Thus the potential modula- 
tion leads for each state with central energy E to side 
bands at energy E + ITiw describing carriers which have 
absorbed I > modulation quanta or have emitted I < 
modulation quanta Tuo. Here we have assumed that all 
potentials oscillate in phase. If one allows for a different 
phase 4> a for each contact a that will add a term e~ tl ^ c " 
to each term in the sum in the wave-function above. Be- 
low, for simplicity, we assume that all contact potentials 
are in phase. 

We now suppose that the modulation potential exists 
only far away from the conductor and that the modu- 
lation potential vanishes as we approach the conductor. 
Thus there is a transition region from a portion of the 
lead in which the potential is oscillating and a portion of 
the lead close to the conductor where we initially assume 
that the potential is time-independent and equal to the 
equilibrium potential. Now we need the wave function in 
the time-independent potential region. This leads to a 
matching problem. If the transition is adiabatic a state 
with energy E in the conductor obtains a contribution 
from all reservoir states with central energy E — ITioj due 
to its side band of amplitude J; (t^O at energy E. In the 
notation of second quantization the annihilation operator 
of an incoming state close to the conductor is 



l (E) = '£k , a!n (E-mu)j l 



hu) 



(13) 



up to corrections of the order of Uui/Ep which arise 
from t he difference o f the wave vectors of the sidebands 
pi = y / 2m(E + lhuj)/h and the wavevector at energy E. 



The current operator Eq. ( |11[ ) is expressed in terms of the 
incoming (and outgoing) states of the stationary time- 
independent scattering problem. Eq. ( |l3"| ) can now be 
used to find the current operator in terms of the reser- 
voir states a.' a n . The current operator becomes 



I a (t) = ~ I dE I dE 



ry$ Ik— — GO 



k')UE-inw) 



A 7 j(a, E, E')a! 5 (E' — khtu) (14) 

where we have introduced the current matrix |?9| 

A Sy (a, E, E') = 5 aS 5 ai l a ~ sl s (E)s a ^E'). (15) 

It is assumed that the modulation imposed on the sys- 
tem is so slow that the contacts can still be regarded as 
being in a dynamic equilibrium state. Thus the quantum 
statistical average can be found by evaluating averages of 
the k a (E — Ihui) as for an equilibrium system. Replacing 
the SL a (E — Ihu)) by their equilibrium statistical expecta- 
tion values we find, 

I a {t) = r d£'^TrA 77 (a, E, E + (k - Z)M x (16) 



■j, Ik 



J, 



eV 7 
Tiuj 



Jk 



eVy 



7 \ -i(k-l)u/t 



f 7 (E-lhuj). 



where f-y(E) = f(E — /j,~ f ) is the Fermi distribution func- 
tion for contact 7. Here /i 7 is the electrochemical po- 
tential of reservoir 7. In Eq. ( |l6[) the trace is over all 
channels in lead a. Taking into account the symmetry 
properties of the current matrix under exchange of the 
energy arguments it can be shown that the current given 
by Eq. ( |l6|) is real. 

From Eq. (|l^) we find that for the dc-current only 
the terms I = k contribute. In this case, as is seen by 
looking at Eq. (fl6|), the energy arguments of the current 
matrix are equal. The trace of the current matrix at 
equal energy arguments and equal lower lead indices are 
just transmission and reflection probabilities. We define 
T aj (E) = —TrAyy(a,E,E). For unequal indices a and 
7 this is the transmission probability for carriers incident 
in lead 7 to be transmitted into contact a. If also a = 7 
the trace of the current matrix is equal to the proba- 
bility R aa of carriers incident in lead a to be reflected 
back into lead a, minus the number of quantum channels 
N a at energy E. In this notation, particle conservation 
in the scattering process is expressed by the sum rule 
J2-y Tay = 0. For the dc-current we find thus 



la(0) 



J dE^T ai {E)jf 

7,1 



eVj_ 

TlLU 



f^E-lluJ). (17) 



Now we expand in this expression the Bessel functions 
in powers of the applied oscillating potentials V 1 . The 
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zero-th order terms gives the dc-current / dc '(°) [{ Vp}] that 
flows as a consequence of stationary differences in the 
applied potentials. We use a superscript (0) to denote 
a response to an external potential only. Since the po- 
tential in the interior is kept fixed this I — V character- 
istic is not gauge invariant. A discussion is provided in 
Rcf. |g9[ and by Christen and one of the authors f35| . 
The next term is second order in the amplitudes of the 
oscillating voltages. For identifying conductance coef- 
ficients recall that the applied potential is of the form 
Vy(t) = ±V 7 {u)e iuJt + iy 7 *(o))e- iwt , with the amplitudes 
taken as real. Thus, from the calculated response to V 1 {t) 
we need to extract the response to the Fourier ampli- 
tudes. These second order terms are determined by the 
photo-conductances 

g^l i [ L o,-u-{V s {Q)}]=-5f]~ < ^ J dET aP [E;{V s (0)}] (18) 
fp{E + M + f (E- M - 2f (E) 

(M 2 

The photo-conductance g„p~[w, — oj; {V£(0)}] determines 
the zero- frequency current in contact a in response to a 
second order voltage oscillation Vg(u>) at contact (3. Note 
that the external photo-conductance generated by bilin- 
ear products Vp(ui)V 1 (uj) with /3 unequal to 7 vanishes. 
Instead of a second order difference in Fermi functions 
we can express the photo-conductance as a second order 
difference of transmission probabilities 

9%\», ~w] = -5 Pi e - J dEf p {E) x (19) 
T a0 (E + hu) + T Q0 (E - Tito) - 2T aP {E) 

For simplicity we have not explicitly indicated the de- 
pendence on the stationary potentials 1^(0). Eq. ( |ig| ) 
shows clearly that we obtain an externally induced photo- 
current only if the transmission probabilities through the 
sample are energy dependent. Thus, for a quantum point 
contact or for a quantized Hall conductor, where we en- 
counter situations characterized by transmission proba- 
bilities which are either zero or one, there is no exter- 
nally induced photo current. This form of the photo- 
conductance also makes it evident that current conser- 
vation is satisfied due to the unitarity of the scattering 
matrix: The sum of all photo-conductances over all con- 
tacts adds up to zero, ^2 a [uj, — w] = 0. However, 
similar to the dc I-V characteristic these conductances 
are not gauge invariant. The sum ffi^I^i — w ] does 
not vanish and consequently the photo-current evaluated 
with these expressions depends not only on voltage dif- 
ferences. 

Let us next consider the displacement current. The 
current at the frequency to is determined by the terms in 
Eq. ( |l6| ) for which k — I — 1 and it is given by 

I a (uj) = ^[ TrA 77 (a, E,E + Tuj) x 

1,1 



Linearising the response to an oscillating external poten- 
tial yields the admittance previously found p9|j30[ |, 

2 p 

5kW)}] = X / dETrA^[a,E,E + fkj;{V 7 (0)}] 

y J {E)-ME + hij) _ 
Tim 

The external admittance given by Eq. ( pl| ) has been the 
starting point of a self-consistent discussion of ac trans- 
port based on the scattering matrix approach. The ap- 
proach has been illustrated in a number of works [[fO] [42| . 
The next term in the expansion is third order in the os- 
cillating potentials and will not be needed here. 

We remark that the external photo-conductances 
Eq. ( |T9| ) are like the dc-current determined by transmis- 
sion probabilities only. In contrast, the displacement cur- 
rent invokes products of scattering matrices at different 
energies and thus depends also on the phases of the scat- 
tering matrix. Expressed in a more physical language, 
the displacement current is sensitive to the densities of 
carriers, expressed here via energy derivatives of phases. 
Below we find that the self-consistent photo-current con- 
tains in fact not only transmission probabilities but, like 
the displacement current, also information on the charge 
accumulated in the conductor. 

Before considering the effect of screening, we discuss 
the relation of the external response to previous work. A 
discussion of shot noise in a conductor with applied ac 
voltages can be found in Appendix A. 



A. Two-terminal conductors 

We consider a two-terminal conductor which consists 
of a tunneling barrier connected on either side to a large 
contact. This is the geometry considered by Tien and 
Gordon @ and Tucker @. The results we obtain from 
the external response described above are in agreement 
with these earlier works. Below we emphasize the arbi- 
trary nature of these results arising from the lack of a 
selfconsistent treatment of the Coulomb interaction. As 
a consequence, different but physically identical configu- 
rations of voltages lead to different results. Later we will 
show how the results for these two arrangements change 
when screening is taken into account, for the specific ex- 
ample where the barrier is a resonant tunneling barrier. 

First we consider the Tien and Gordon case, where 
one of the contact potentials is oscillating and the other 
is kept fixed, Vi(co) = V(uj) and V-z{<jS) = 0. For sim- 
plicity we assume that the scattering matrix has been 
diagonalized such that transmission through the barrier 
is described by a transmission probability T m {E) and a 
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reflection probability R m (E) for the m-th eigen channel. 
Using Eq. jl7|) and using the sum rule for Bessel func- 
tions, J2i Ji+k{x)Ji{x) = S k o, we find 



MO) 



-IjdE^Jf 



eV(u) 

TluJ 



T m (E) x 



x [h{E + ZM-/2 (E)}. 



(22) 



Tien and Gordon § express the transmission prob- 
abilities with the help of Bardeen's formula T — 
An\t\ 2 i'i(E)i'2(E + ITiLo) in terms of a matrix element t 
and the density of states (E) and v 2 (E) to the left and 
right of the barrier. In our work the energy E is a global 
variable, whereas Tien and Gordon measure energy in 
the densities of states from the conduction band bottom 
to the left and right of the barrier. 

The time-dependent current was investigated by 
Tucker for the same geometry 0]- Using Eq. ( |l6| ) we 
find 

J x (i) = ^ TrA 77 (a, E,E + khuj) 



fko 



xf 7 {E- Ihu). 



Too 



(23) 



This is the result of Tucker Q except that he considers 
the barrier transmission to be energy independent. 

However, the results obtained above are not invari- 
ant under an equal shift of all potentials. For example, 
in the Tien and Gordon case, an experimentally equiv- 
alent situation would be to set Vi(uj) — V(oj)/2 and 
V2(o>) = —V(u>)/2. This, however, yields a different 
result in the non-interacting theory. Even worse, set- 
ting Vi(tu) = V2(u>) = V(ui)/2 should yield no photo- 
current, but gives the same as for Vx(u>) = V(u>)/2 and 
= — V(ui)/2. To remedy this we introduce in the 
next section a simple selfconsistent scheme to achieve 
charge and current conservation, similar to one used pre- 
viously mm. 



N= J2 fd 2 r[ dE^(E)^(E) x 

af3nm 

K„ (r, E)y pm (r, E)al n {E)a 0m (E) , 



(24) 



where v an {E) is the density-of-states for channel n in 
contact a, and ^ an {r,E) is the corresponding wave- 
function for a scattering state describing carriers incident 
in contact a in channel n. 

We now define the partial density-of-states matrix 



dE 



with elements 



ijE = J dh^(E)^(E)K n (r,E)y 0m (r,E). 

(25) 

This matrix can also be expressed in terms the scattering 
matrix and its derivatives IH3] 



dM 1 s 
dE 



TT1 t-^ 



dE 



dE 



sps(E) 



(26) 



Using this and Eq. ( |l3[ ) we find the number operator in 
the presence of oscillating contact potentials 



N=J2jdEJ l (^)J l 



afjlk ' 



huj 



(k'tiE-lh^^k^E-khu;), 



dE 



with the expectation value 

'eV a 



N 



J^JdEJf 



(27) 



(28) 



We can shift the frequency dependence from the Fermi 
function to the partial density of states 



dE 



(29) 



B. Density operator 

When applying voltages to the conductor, the sample 
will be charged. The net charge of the sample in response 
to a potential applied to a contact can be decomposed 
into two contributions: A charge response, called the in- 
jectancc of the contact, at fixed internal electric poten- 
tial and a charge response due to an electrically induced 
potential. Here we determined the injectances of a multi- 
terminal conductor. In the next section these results are 
used when treating the problem with interactions. 

At zero frequency, the number of electrons in the sam- 
ple is determined by the operator Mm 



dN (0) 

and thus identify the injectance — jg— at energy E in the 
presence of a potential variation at contact a. Here the 
upper index is once more used to emphasize that this 
density is evaluated at fixed internal potential. To second 
order in the oscillating potential, V a (u>), the injectance is 



dN^ 



dE 



Tr 



dM aa (E) e 2 |U Q H| 2 



dE 2 {Tiuj) 2 

f dN aa {E + Huj) ^ d\f aa (E - huo) 



(30) 



dE 



dE 



, cW aa (£) 
dE 



In the limit that |V a (u;)| becomes small compared to Tiio 
the injectance is that produced by a static voltage. 
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III. INTERNAL RESPONSE: SELF-CONSISTENT 
SCREENING 

In response to a potential variation at a contact the 
charge distribution in the interior of the sample is driven 
away from its equilibrium pattern. Coulomb interactions 
oppose such a variation. In the problem of interest here 
a variation of the sample charge can come about both 
because we in general consider a biased sample such that 
a dc-current flows and because we subject the sample to 
ac-voltages. In general it is a non-equilibrium dynamical 
potential landscape that matters. Here for simplicity we 
consider the sample to be zero-dimensional and assume 
that it suffices to consider a single internal potential U. 
Such an approximation is often used in the literature on 
the Coulomb blockade and in the scattering approach to 
electrical conduction has been used to discuss the non- 
linear I-V characteristic of mesoscopic samples |3q| and 
ac-transport in Refs. [ ^3|j3C|j45[ | . At equilibrium, if all 
voltages at the contact of the sample are equal, and in 
the absence of ac-potentials, the value of this potential 
is U = U eq . Our first task is to determine the zero- 
frequency part of this potential. 

To be more specific we now consider a sample coupled 
to a gate, as an example see Fig. 1. We denote the con- 
tact to the gate by the index g and the capacitance of the 
central region of the conductor to the gate by C g . The 
capacitance between the central region of the conductor 
to the reservoir a is denoted by C a . Next we introduce 
an index v which runs over all N current contacts of the 
sample v = a — 1,2, ..N and in addition includes the 
contact to the gate v = N + 1 = g. 



A. Static internal potential 

Consider first the equilibrium potential t/g 9 . The 
grandcanonical potential with the Coulomb energy in- 
cluded is minimal for a potential V^ q that obeys the 
Poisson equation. In our case the Poisson equation is 
discretisized and is expressed with the help of the geomet- 
rical capacitances introduced above. The net electronic 
charge on the sample is that permitted by the Coulomb 
interaction: 



Q-Q + = ^C V {U-V V ) 



(31) 



Here Q is the electronic charge, Q + is an effective "ionic 
charge" created by the donors and C„ are the geometrical 
capacitances. 

For V v = 0, the equilibrium charge Q = Qq 9 and the 
equilibrium potential U — follow from Eq. (|3l]) as 
follows. The electronic charge on the conductor can be 
expressed as a sum of all the charges injected from the 
various contacts, 



r q 



E 



where the injectance [^5|33| of contact a is given by 
Eq. (30). Note that the scattering matrix and thus the 
injectance also depends on U^. Eq. ( |32] ) is thus a self- 
consistent equation for the equilibrium potential. 

Next, let us keep the ac-voltages turned of but apply dc 
voltages to the contacts, charge will flow into the conduc- 
tor causing a shift of the static potential in the barrier. 
We denote the resulting potential by Uq. It is a function 
of the applied potentials V v since now the injected charge 
depends on all the applied voltages. The injected charge 
is given by 



Qo = E/ 



dE 



(33) 



Using Eq. (|3^) to express the effective background charge 
in terms of the scattering matrix and the charges on the 
capacitors gives the following self-consistent equation for 
determining the internal static potential in the sample 



E / dE 

= ^, c A u o - v v ). 



E 



/_Njm dE 

dE 



(34) 



(32) 



Here V a = fi a — /J, is the deviation of the electrochemi- 
cal potential in contact a from its equilibrium value \x. 
This approach was used by Christen and Biittiker [[35| to 
study the nonlinear conductance for a resonant tunneling 
barrier. 

Next, consider next the case that is really of interest 
here. In addition to possible static voltage differences 
we have time-dependent potentials at the contacts. As a 
consequence the (unscreened) charge Q(t) in the sample 
is also a function of time. It can be Fourier transformed, 
and we expect Fourier components at the oscillation fre- 
quency uj of the voltage and at all higher harmonics, kui. 
As a consequence the potential inside the conductor will 
also oscillate and will similarly have Fourier components 
at all harmonics, U(kuj). If an oscillating voltage at a 
contact, due to non-linear processes, also changes the 
time-averaged charge in the sample then the potential Uq 
as determined above would be modified by the presence 
of the oscillating potentials. To take this into account we 
write the injected charge as a response to external po- 
tentials in the presence of a self-consistently determined 
static potential plus the response from the internal oscil- 
lating potential. The response to the internal potential is 
determined by three unknown response coefficients, Xiai 
Xai and Xu such that, 

/dN^ f 
e-^-dE + J^ / dE Xia (E)U*(<j)V a (cj) 

+ J2 J dE Xa i{E)V*{u)U{u)+ X ii{E)\U(u)?. (35) 
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To determine \ia, Xai and Xu we use the fact that the 
injectance be invariant under a shift of all oscillating po- 
tentials by an equal amount. This yields the coefficients 



Xia(E) = X<*i(E) = Tr 



1 r_e_\ 2 ( dU aa {E + Tiuj) 
2\h^J V dE 



dE 

X»(E) = -J2xUE). 



dN aa {E - hu) _ 2 dN aa {E) 



dE 



(36) 
(37) 



With this we can express the gauge-invariant injectance 
as 



dN a 
~dE 



Tr 



dM aa (E) 
dE 



\V a {u)~U(uj)f 
(HUJ) 2 



(38) 



dM aa {E + huj) : dM aa {E~Tiuj) 2 <lN aa {E) 



dE 



dE 



dE 



Eqs. (34) and ( p8| ) now allows us to find the static 
internal potential Uo in the presence of static and os- 
cillating contact voltages. Note that ^jg- depends on Uq 
since the scattering matrix depends on Uq. The potential 
C/o depends on the dc-voltages applied to the sample and 
depends through non-linear processes on the amplitudes 
of the ac-voltages and the frequency. 

Our next task is now to find the current response to 
the oscillating internal potential U (t) . 



B. dc-current 



#TO(o)}] = WW}] WW 

01 

+ ]T g aPi [u, -to; {V B (0)}]Vp(u)U*{u) 



+ 53fl«frk -w; {V s (fl)}]U{u)V;{u) 

7 

+g aii [cj,-u;{V s (0)}]U(uj)U*(ui) (40) 

where the index i refers to responses due to the internal 
potential U(u>). 

The responses to the internal potential are found by 
demanding that the current is invariant with respect to 
a shift of all voltages (gauge invariance). Lowering all 
voltages at frequency u> by U(u>) shifts the internal poten- 
tial to the external potentials. Comparing the resulting 
expression with Eq. ( ^p| ) gives 

Wk-«;W)}] =-£0$>,-w;{V f (o)}], (41) 

7 

gc*r[u,-u;{Vt(p)}] = ~E {V s (0)}], (42) 

P 

g aU [u,-u J ;{V s (0)}} = ^g i Xl",-u;{V S (0)}}. (43) 
0-1 

Thus the photo-responses to the internal potential 
are determined by combinations of external photo- 
conductances. 

With these conductances, the dc-current can be writ- 
ten as 

/«(0)=l£TO(0)}] (44) 
+ Y,9%["-^{V i m\\V0{u)-U{uj)\ 2 . 



Consider first the photo- induced dc-current. The dc- 
current can be divided into two parts, one due to direct 
transmission processes, and one due to transmission after 
absorption (emission) of a photon followed by its emis- 
sion (absorption). Both processes take place in a self- 
consistently determined electrostatic background, which 
depends on all voltages at all frequencies. 



i a (o)=i* c [{v (o)}] + it [Wo)}] 



(39) 



Here ^ c [{V/3(0)}] is determined from the first term of 
the sum in Eq. (|l7|), where now the scattering matrix 
depends on Uq. 

The photo-current can be written generally as the sum 
of the response to the external oscillating potential and 
the internal potential U(u>). To proceed we now consider 
e = eU (u>) I (hu>) , a small parameter in which we can ex- 
pand. All the oscillating contact potentials are also of 
order e. In this work we will stop this expansion at the 
first non-trivial order. Since photon-assisted tunneling is 
of second or higher order in the oscillating potentials, we 
carry the expansion to second order. 

For the photo-current we obtain 



Note that the current depends on the difference between 
an applied voltage and the internal voltage only. 

All the non-linear transport coefficients in Eqs. ( ^0[j4^ ) 
also depend on Uq, the self-consistent dc-potential. 



C. Displacement current 

The current at frequency u> is only needed to first or- 
der in the applied oscillating voltages. In addition to the 
external potential the oscillating internal potential also 
contributes to the current. In the presence of the internal 
potential the general form for the current is to first order 
in the potentials 

,(°)r 



W=£ffi>;{W}]W 



+g ai [L,;{V 1 (0)}]U(u). 



(45) 



Here <7^a[w; {V^(O)}] are the external ac-conductances 



given by Eq. ( J21| ) and g ai is the ac-response to the in- 
ternal potential. Again we determine g a i through the 
requirement that this expression is invariant under an 
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overall shift of the potential. This gauge invariance ar- 
gument determines the response to the internal poten- 
tial in terms of external responses; g a i[u>', {V 7 (0)}] = 

-E^SkWo)}]. 

Both for the de-current and the ac-current we now 
know the response to the external voltages V 7 (0) and 
to the internal potential U(t). But the internal potential 
is thus far not determined. This is our next task. 

To be more specific we now return to the sample shown 
in Fig. [j]. The current at contact a is the particle cur- 
rent plus the displacement current (capacitive) current 
I a (uj)— iuC a (V a (ijj) — U(u>)) with I a as determined above. 
The current from the gate to the sample is purely capac- 
itive and is given by I g (u>) — —iuj(V g {uj) — U(u>)). Since 
the overall charge at frequency u> is conserved the sum of 
these currents must vanish. Thus we must have 



^j tt H = - 1 ^c„([/H-v;H). 



(46) 



Solving this equation for the internal potential yields 



E Q/3 sSM^(o)}]-^£^ 



(47) 



The external ac-conductances and the geometrical capac- 
itances determine the potential U(uS) and determine the 
self-consistent dc-current due to photo-assisted tunneling 
and the self-consistent ac-conductances. 



IV. RESONANT TUNNELING BARRIER 

As an application of the selfconsistent theory devel- 
oped above, we consider the photo-induced dc-current 
through a resonant tunneling barrier. The experimental 
setup is taken as sketched in Fig. |l|. Each side of the 
resonant barrier is connected to reservoirs with chemical 
potentials /ii and /i2 and capacitances C\ and C%. The 
interior of the barrier is coupled to a gate with a capac- 
itance C g . For simplicity we assume that the gates are 
macroscopic with no dynamics of their own. A dc bias 
will be applied by making eV = e(Vi — V2) = Mi — ^2 
non-zero. 

The scattering matrix close to a resonance is given by 
the Breit-Wigner formula 



E-Ea- eUn 



*7/2. 



(48) 



Here 7 „, n = 1, 2 are the partial widths of the resonance 
proportional to the tunneling probability through the left 
and right barrier and 7 = J2 n 7« ^ s the total width of 
the resonance. 5 m are the phases acquired in the reflec- 
tion or transmission process and Eq is the position of the 
resonance. The term eUo = e ¥lM±¥zi21 + \y ensures 



invariance upon a shift of the dc voltages p5fl. W is de- 
termined by the condition Eq. ([54]) and Eq. (p8|), and is 
a function of V~i(0) — V^O) only. The injectivities are Ej| 



1 



la 



dE 2tt (E-E - eU ) 2 + ( 7 /2) : 



(49) 



The Breit-Wigner formula is a reasonable form for the 
scattering matrix as long as the energy does not get close 
to the next resonance level. Assuming that the level spac- 
ing of our system is large enough such that neighboring 
levels can safely be ignored we will use the formula in a 
wide energy range. 




FIG. 2. Differential conductance as a function of dc bias 
from the non-interacting discussion. The left contact poten- 
tial is oscillating. The parameters are hu>/(j/2) = 5, e 2 = 0.1, 
and the Fermi energy is equal to the resonant energy /j, = Eo . 
For comparison the dashed line shows the transmission prob- 
ability determined from the Breit-Wigner expression. 

Photon-assisted tunneling is most easily seen either in 
the differential conductance as function of bias voltage 
dI(V)/dV, where side-peaks shows up at multiples of 
the photon energy, or in the dc-current when varying the 
gate potential [gf I(V g ). In Fig. ||we show an example of 
a dI(V) I dV-curve using the non-interacting discussion, 
Eq. ( p2| ) and using the Breit-Wigner expression Eq. (|i^ ) 
with Uq = 0. The potential of the left contact oscillates. 
We apply a dc voltage V = VI — V2, take hui/(j/2) = 5 
and consider the symmetric case 71 = 72 = 7/2. In 
this and all the following examples we use e 2 = 0.1, for 
which the expansion to second order is pertinent. For 
this choice of parameters only the first side band peaks 
can be resolved. In a non-interacting discussion one iden- 
tifies Uq — V g and as a consequence both the differential 
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conductance as a function of dc-voltage, or the current 
as a function of gate voltage, the side bands are observed 
at a voltage corresponding to the photon energy. Dis- 
cussions which neglect interactions do not discriminate 
between these two methods for observing photo-assisted 
transport. If we now consider the physically meaningful 
result, the theory which includes interactions, the gen- 
eral behavior will remain the same, but the two methods 
of analyzing photo-assisted transport, i. e. considering 
dI(V)/dV or I(V g ) now give in general different results. 
The effects brought about by screening, discussed in more 
detail in the next section, are: First, the relative weight 
of the sidebands and the central peaks will not be the 
same in the two situations. Second screening also brings 
about an asymmetry in the weights of the side bands 
for ±nhuj. In a discussion that neglects interaction the 
side bands have the same weight. In contrast, in the in- 
teracting case, if the equilibrium chemical potential does 
not coincide with the resonant energy, screening will be 
different for the two voltages where peaks are seen, and 
accordingly their weights will differ. Such asymmetries 
are seen in experiments |7j . Below we discuss these effects 
in detail. 



A. Gate driven case 

First, consider a sample subject to a dc bias V — 
V\ — Vi and an oscillating voltage V g (u>) applied solely 
at the gate. For simplicity we take C\ = C2 = 0. In this 
case there can be no dc photo-current when fix = /i2, 
since J2p A /3/3( a > E, E) 



(see Eqs. QL8|) and (|44|)) as 
a consequence of the unitarity of the scattering matrix. 
The effect of photon-assisted tunneling in this setup is 
controlled by the internal potential. Thus, it is of in- 
terest to understand how it relates to the applied gate 
voltage in the presence of screening. From the selfcon- 
sistent theory (see Eq. ([!?])) we find for the ratio of the 
applied to the external potential 



Vg(u>) 



1 £*>'*0 

a/3 



jC 



(50) 



This ratio is determined by the ac-conductanccs 
gj^p(u>;V). These ac-conductances are known. At zero 
temperature, for the symmetric resonant tunneling bar- 
rier 71 = 72 they are given by Fu and Dudley p9| and for 
the asymmetric case 71 7^ 72 by Biittiker and Christen 
il: 



(°)/ >l (°){ \ 

9n M =32i M 

(0)/ \ (0) / A 

322 M =012 H 



gi ( ?(u;V)=g ( 3(w;-V), 
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7i 
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(51) 

(52) 
(53) 



3i° 2 ) M 



e 2 7172 



h 7?kj 1 



(54) 
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i [(ft -hm-E -W~ eV/2) 2 + ( 7 /2) 2 ] 
2 n [fa-Eo-W-eV/2)* + {>Y/2)*] 

1 [(p + huj - E - W - eV/2) 2 + ( 7 /2) 2 ] 

2 n [(/x - E - W - eV/2Y + ( 7 /2) 2 ] 



arctan 



arctan 



H + huj - E - W - eV/2 
7?2 

H - fiu - E - W - eV/2 
7/2 



With these expressions Eq. (jMj), the ratio of internal 
to external potential, can be evaluated. This ratio has 
two simple limits. In the non-interacting limit C — * 00, 
the internal potential directly follows the applied poten- 
tial. In the limit C — > 0, we have a charge neutral sample 
and U(lj) = 0. 
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FIG. 3. Ratio of the internal potential to the gate volt- 

2 

age as function of the Fermi energy, for C — ^ , V = and 
for the frequencies a) — 10, 6) fki>/(pj/2) — 3, and 

c) Rw/( 7 /2) = 1. 

In Fig. ||, we show the absolute square ratio of the in- 
ternal to the external potential for different frequencies, 
when sweeping the Fermi level through the resonance. 
The non-screened case C — > 00, where the ratio is 1, is 
shown as the dashed line. It is evident that screening in- 
troduces a large renormalisation of the internal potential 
for this choice of capacitance with a strong dependence 
on frequency. One observes the largest effect when the 
Fermi energy is close to the resonance. This is expected 
since the density in the barrier is a Lorentzian with a 
peak at resonance [M, thus providing more screening 
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electrons. As a function of frequency the ratio changes 
qualitatively; for low frequencies the internal potential 
is reduced compared to the external potential, whereas 
with increasing frequency the situation reverses. 
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FIG. 4. Ratio of the sideband weight to central peak 
weight as function of capacitance in the current versus gate 
voltage characteristic I(V g ) for frequencies a) hw/( , y/2) = 3, 
b) &d/(7/2) = 5, and c) Tilu/^/2) = 10, when e = 0.1. The 
dashed line shows the result when no screening is present. 
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FIG. 5. Ratio of the sideband weight to the central peak 
weight as function of the the capacitance in the differen- 
tial conductance dl/dV for the frequencies a) hui/( , y/2) — 3, 
b) hw/(j/2) = 5, and c) hu>/(-y/2) = 10 and e 2 = 0.1. The 
dashed line shows the result when no screening is present. 



Next consider a the current as a function of gate volt- 
age. Since screening depends on the position of the reso- 
nant level compared to the equilibrium electro-chemical 
potential, the central peak and the sideband will expe- 
rience a different degree of screening and, thus, their 
weights will no longer be given by a Bessel function be- 
haviour as in the non-interacting theory. In Fig. ^ the 
ratio of the sideband peak to the central peak is shown. 
The non-interacting theory predicts a ratio of 0.125 for 
the parameters chosen (dashed line). It is seen that, de- 
pending on capacitance and frequency, this ratio can be 
quite different. 

Similarly, when measuring the differential conductance 
as function of dc voltage screening will also vary as func- 
tion of voltage. In this case the sideband weight to central 
peak weight ratio is shown in Fig. Again, large differ- 
ences with respect to the non-interacting case is possible. 

An interesting effect due to the dependence of screen- 
ing on the dc voltage (or the gate voltage) is that side- 
bands will no longer be strictly Lorentzian, but skewed. 
However, this skewing effect is rather small and probably 
difficult to resolve experimentally. 



When the Fermi level is off resonance the first side- 
bands corresponding to absorbing and afterwards emit- 
ting a photon and visa versa occur at different poten- 
tials. Screening will therefore occur asymmetrically for 
the two peaks introducing an asymmetry between the ± 
sidebands. This effect is illustrated in Fig. [| Experi- 
mental observation of this effect has already been made 
jjj , although it has not been studied systematically. 



Another effect is visible in the inset in Fig. |[ One 
notices that the width of the central peak is significantly 
larger than the width of the sidebands. Since the ca- 
pacitance in this example is rather small, the charging 
energy is large, and when increasing the dc bias voltage 
the added charge gives rise to a huge increase in the static 
internal potential. The result is that the resonance floats 
upwards in energy, widening the peak. For the same rea- 
son, the distance from the central peak to the sideband 
is no longer simply hu>, but substantially larger. 
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FIG. 6. Weight asymmetry for the ±7kj sidebands 
as function of capacitance in the dl/dV-characteristic 
for 71/(7/2) = 1/4 and 72/(7/2) = 3/4, and for 
a) hw/(j/2) = 3, n/(rf/2) = -5, and 6) huj/(j/2) = 5, 
/i/(7/2) = —10. The dashed line shows the result when no 
screening is present. The insert shows an example of a differ- 
ential conductance curve as a function of the bias voltage for 
Tl = 72 , Rw/( 7 /2) = 5, n/(rt/2) = -10, and C = 0.l£. The 
dashed line is the result without screening. 



B. Contact driven case 



Since g l2 (oj) / (ji^) is a function only of 7 we find that 
the zero bias current is proportional to the effective 
asymmetry of the double barrier, 71 ~ 72 . For small ca- 
pacitances, Cji^2/(^9i2 ( w )) ^ 7ii72) t ne current is 
directly proportional to the asymmetry of the barrier 
without any renormalisation from screening Ii(Q;V — 
0) = gm(uj,-u)\V 2 (uj) - Vi(^)| 2 ^i. Thus, the ac- 
field effectively pumps electrons through the system. The 
noninteracting Tien and Gordon discussion, as given by 
Eq. (p2|), also predicts a current at zero bias. The re- 
sult, however, being independent of the asymmetry of 
the system. This prediction of a zero-bias current for 
the symmetric case is due to the violations of gauge- 
invariance. That a symmetric structure, in the absence 
of dc- voltages, cannot exhibit a photo-current, can be 
understood from the following symmetry and invariance 
conditions. Consider first a variation of the voltage at 
the left contact V\(lo) = Vb cos(wi) and suppose this pro- 
duces a dc-photo current 1%. Then consider a voltage 
variation of the right contact Va(w) = — Vbcos(urf) = 
Vb cos(ujt + 7t). By symmetry this must give a current 
I2 = In reality, however, due to gauge invariance 

these two voltage oscillations are experimentally the same 
and hence must give rise to the same dc-current. But the 
only dc-current which reflects this symmetry is I = 0. 
Clearly, the correct answer is a consequence of gauge in- 
variance. 



A setup often used experimentally is to couple the os- 
cillating field to the conductor via a bowtie antenna 0. 
In this case we assume that there is no gate, C g — 0. For 
simplicity we take the capacitances across each tunneling 
barrier to be identical, C\ = C2 = C/2. The dc-current 
into contact 1 is then 



gfh{u,-u);V) 

+9 i 122( UJ 1 -UJ]V) 

|^ 2 H-Fi(w)| 2 



(55) 



For the symmetric case 71 = 72 the dc current will van- 
ish for zero bias, because of the symmetry of the prob- 
lem. However, in contrast to the gate driven case a zero- 
bias current can be generated for the asymmetric sample, 
given by 

h(0; V = 0)= flSl(w,-o;)|y 2 (w) - VI H| 2 ^^ 

7 



V. CONCLUSION 

We have extended the scattering matrix approach to 
transport in phase-coherent conductors to take into ac- 
count oscillating contact potentials and internal poten- 
tials in nonlinear order. The effect of screening has been 
taken into account to second order in the oscillating po- 
tentials by means of an RPA treatment. The result is 
a theory, valid for arbitrary dc voltages, which is cur- 
rent and charge conserving (gauge invariant). The in- 
ternal potential in the conductor has been treated as a 
single parameter. Certainly, to go beyond this approx- 
imation and treat a more realistic continuous potential 
distribution would be interesting. But even for the case 
of linear ac-transport, a scattering matrix for continuous 
potentials exists only to linear order in frequency |jl"|| , 
and exceptionally to second order Q). Discussions of 
the dynamic conductance of a ballistic wire over a wide 
range of frequencies, taking into account spatial potential 
variations |32|, are not yet formulated within the scat- 
tering approach. It would also be interesting to extend 
our discussion to higher order in the applied voltages. 
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For large field strengths it is possible to make one of the 
Bessel functions zero, giving rising to dynamic localiza- 
tion |5(]||[l2|]. Since we find that Bessel functions can in 
general not give a gauge invariant answer it is clear that 
the criteria for dynamic localization will be changed in 
an essential way in the presence of interactions. 

We have applied our theory to photon-assisted tun- 
neling using a resonant tunneling barrier as an example. 
The two standard setups for photon-assisted tunneling, 
applying the modulation to one of the contacts in a two- 
terminal experiment, or coupling the potential to the con- 
ductor via a gate was examined within the self-consistent 
theory. In both cases, the inclusion of screening leads 
to a renormalisation of the non- interacting answer. The 
driving field is not the applied field but the total field. 
Since the effective field is dependent on screening and 
therefore on applied bias, chemical potential etc., the 
weights of the central peak and the side peaks in the 
differential conduction versus applied voltage differ from 
the non-interacting theory. Furthermore, the peak weight 
is no longer distributed according the increasing order of 
Bessel functions. This leads to the peak ratios being a 
complicated function of the screening properties of the 
system, and predicts an asymmetry between the corre- 
sponding left and right sidebands. Asymmetric photo- 
conductance peaks have been observed |?]||. 

The necessity to include screening in the treatment of 
photo-assisted transport is most clearly exemplified by 
the following consideration. For a spatially symmetric 
system a non-interacting theory (Tien and Gordon) pre- 
dicts a photo-current in response to the oscillation of ei- 
ther the left or the right contact voltage. In contrast, the 
gauge invariant discussion presented here, predicts that 
a symmetrical system exhibits no photo-current. Our 
result for the two-terminal resonant tunneling barrier, 
Eq. (|5^) is a photo-current which is proportional to the 
asymmetry of the tunneling rates of the resonant double 
barrier structure. 



In this work we have emphasized that interaction ef- 
fects are important whenever a variation of a parameter, 
an oscillation of a voltage, changes the charge away from 
its equilibrium value. In photo-assisted tunneling it is not 
sufficient to consider just the dc-current, but a theoreti- 
cal discussion has to be self-consistent at all frequencies. 
Thus there is necessarily a relation between the photo- 
assisted dc-current and the displacement current. Only if 
the charge is investigated at all frequencies can an elec- 
trically meaningful, that is gauge invariant, answer be 
found. 
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APPENDIX A: CURRENT NOISE 

The analysis of this paper concentrates on the average 
zero-frequency photo-current. However, the approach 
used here also allows to find the fluctuations of the cur- 
rent. Of particular interest are the current-current cor- 
relations which determine the spectral densities of the 
current fluctuations. Here we present the general result 
for the noise spectra of a multi-terminal conductor in the 
presence of oscillating contact potentials assuming that 
the internal potential is kept fixed. As with the aver- 
age dc current a physically meaningful result requires in 
general a discussion of the effects of screening. 

For a multi-probe conductor with potentials V a cosfuA) 
at frequency uj applied to the contacts, using Eq. 
we find the correlation function 



Tr[A 7(5 (a, E, E')A Sj (/3, E' + (I' - k)fwj, E + (&' - l)huo)\ 

[/ 7 (E - lhcu)(l - f 5 (E' - khcu)) + f 6 (E' - khio){\ - f 7 (E - Ihu))]. (1) 



Here the brackets {, } denote the anti-commutator. In 
the presence of ac voltages the current-correlation func- 
tion is not only a function of the relative time r but de- 
pends also on the absolute time t. Experimentally what 
is of interest is the noise spectrum on a time scale long 
compared to 2ir/ui. Therefore we define the noise spec- 
trum as an average 



Sap{T) = hj dt{{^L{t + T),Mp{t)}), (2) 

where T = 2tt/lu is the period. The factor 1/2 
arises because we have symmetrized the correlation func- 
tion. The spectral density is related to the current- 
current correlation function via 27rS' , / 3(r2; (v)5(Tl + f2') = 
(1/2)({A/ Q (fi), Al^fJ')}), which is just the Fourier 
transform of S(t). We find 
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- (1)7 « £/< (§) * (§) (§) * (t 

Tr[A 75 (a, E,E + hn)A Sy (/3, E + Titt+{k' - l)hcj, E + (fc' - i)^a;)] 

[/ 7 (£ - lfwj){\ - f s (E + HQ,- khu)) + f s (E + khu)(l - f 7 (E - Ifiuo))]. (3) 

In the limit of vanishing driving frequency, u> = 0, Eq. (|^) reduces to the frequency-dependent noise spectra of Ref. [^9[ . 

For the special case that the scattering matrices can be taken to be independent of energy, i.e. A^g(a, E, E + Tluj) = 
A 7 s{a) Eq- @ simplifies considerably. Using the addition theorem for Bessel functions we find 



eiVs - Vy) 



j6,l 

[/ 7 (E + ITwj){\ - fs(E)) + f s (E)(l - / 7 (E + ITwj))}. (4) 



For a two-terminal conductor this result is identical to 
that of Lesovik and Levitov [E5j even though in that work 
this result was derived in response to an electric field and 
not as here as a response to an oscillating contact volt- 
age. In the experiment of Schoelkopf et al. [2(J the shot 
noise is measured in the presence of an oscillating voltage 
applied to the contacts of the sample. 

We emphasize that the noise spectra given by Eqs. (j^) 
and (^) give only the noise for fixed internal potential. 
We have already remarked that the average dc current 
exhibits an external response due to photo-assisted trans- 
port only if the transmission probabilities exhibit an en- 
ergy dependence (see Eq. (pToj) ) . In contrast, in the shot- 
noise spectra, we have an effect even if the scattering 
matrix is taken to be energy independent. That is a con- 
sequence of the fact that the noise spectra depend in a 
non-linear way on the Fermi functions. 
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